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ABSTRACT 


Exact  methods  have  been  developed  for  calculating  light  scattering  from  particles  on  a  perfectly 
reflecting  surface.  The  computed  solutions  obey  Maxwell's  equations  and  the  appropriate 
boundary  conditions.  A  brief  discussion  of  the  theory  is  followed  by  a  presentation  and  discus¬ 
sion  of  selected  scattering  results.  Comparisons  are  made  with  the  predictions  of  Young's  forward 
scatter  Mie  theory  model.  Under  certain  conditions  the  exact  results  differ  significantly  from  Mie 
theory.  The  exact  differential  aoss  sections  exhibit  large  amplitude  oscillations  that  are  not  pre¬ 
dicted  by  Mie  theory,  llie  bidirectional  reflectance  distribution  function  (BRDF)  is  also  calcu¬ 
lated  using  exact  scattering  theory  for  both  spherical  and  hemispherical  particles  that  are  both 
absorptive  (complex  index  of  refraction)  and  nonabsorptive  (real  index  of  refraction).  These 
exact  BRDF  results  are  compared  to  the  predictions  of  Young's  model. 
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1.  INTRODUCTION 


Light  scattering  from  particle  contaminants  on  the  surfaces  of  mirrors  Is  one  of  the  main  con¬ 
tributors  to  stray  radiation  In  satelUte-bome  surveillance  systems.  This  limits  the  ability  of  the 
system  to  detect  and  track  dim  targets  that  are  near  a  bright  radiation  source  such  as  the  Earth. 
The  designers  of  these  systems  need  reliable  methods  for  quantltadvely  estimating  the  Intensity  of 
the  scattered  light  In  order  to  minimize  the  effects  of  stray  radiation. 

Previous  laboratory  Investigations  of  this  problem  Include  both  experimental  measurements  and 
theoretical  calculations.  The  model  systems  In  these  studies  were  spherical  particles  on  a  con¬ 
ducting  plane.  The  prior  theoretical  calculations  have  been  based  on  approximation  methods  that 
neglect  the  Interaction  between  the  particle  and  the  reflecting  surface. 

Recently,  we  have  developed  and  Implemented  several  exact  procedures  for  calculating  light 
scattering  by  spherical  particles  and  also  certain  types  of  nonspherlcal  particles  on  a  perfectly 
conducting  plane  surface.  It  is  the  purpose  of  this  report  to  present  some  selected  results  of  these 
exact  calculations  and  to  also  compare  these  results  with  the  predictions  of  Young's  forward- 
scatter  Mle  theory  method.^ 
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2.  EXACT  METHODS 


2.1  METHOD  OF  IMAGES 

The  scattering  problem  is  illustrated  In  Fig  la.  A  particle  rests  on  a  perfect  reflecting  (infinite 
conductivity)  plane  surface  mirror.  An  electromagnetic  plane  wave,  represented  by  the  electric 
field  vector  Einc.  is  incident  upon  this  system.  This  gives  rise  to  a  reflected  plane  wave.  Efef.  and  a 
radially  outgoing  scattered  wave.  Escat*  In  the  region  above  the  mirror,  the  field  must  satisfy 
Maxwell's  field  equations  and  also  satisfy  the  condition  that  there  can  be  no  tangential  component 
of  the  electric  field  at  the  reflecting  surface.  The  problem,  as  presently  formulated,  is  quite  diffi¬ 
cult  to  solve  because  of  the  boundary  condition  thtt  must  be  satisfied  at  the  surface.  The  method 
of  images  can  be  used  to  reformulate  the  problem  in  a  much  simpler  form. 


<a) 


Hgure  1.  Illustration  of  the  method  of  images:  (a)  physical  system;  (b)  equiva¬ 
lent  method-of-images  system. 
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An  observer  looking  down  on  the  system  in  Fig  la  will  see  the  particle,  the  mirror  image  of  the 
particle,  the  source  of  the  incident  wave  (e.g.,  a  laser),  and  the  mirror  image  of  the  source.  The 
image  of  the  source  is  the  apparent  origin  of  the  reflected  wave.  This  apparent  picture  is  turned 
into  reality  by  creating  the  system  shown  in  Fig.  lb  where  the  mirror  has  been  removed,  the 
image  of  the  particle  has  been  rqrlaced  by  a  real  particle,  and  the  image  of  the  source  has  been 
replaced  by  a  real  source.  Hie  image  source  is  the  origin  of  the  reflected  wave,  Eref.  which  is 
slmwn  in  the  figure  as  originating  below  the  surface  of  the  mirror  and  incident  on  the  two-particle 
system  along  with  the  original  incident  wave,  Einc-  It  can  be  shown  that  the  solution  to  the  scatter¬ 
ing  problem  shown  in  Fig  lb,  in  the  region  above  the  mirror,  is  identical  to  the  solution  to  the 
pr(rt)lem  shown  in  Fig. la  in  this  same  region.  The  system  shown  in  Fig.  lb  is  easier  to  solve  than 
the  original  system  because  no  explicit  surface  boundar>’  conditions  are  imposed  on  the  solution. 
The  bouiKlary  conditions  at  the  surface  of  the  mirror  are  automatically  satisfied  as  a  natural  con¬ 
sequence  of  the  symmetry  of  the  system. 

SPHERICAL  PARTICLES 

The  problem  is  to  calculate  the  scattering  from  the  two-particle  system  shown  in  Fig.  lb.  If  the 
particles  are  si^ierical,  the  multipole  expansion  method  of  Bruning  and  Lo^  can  be  directly 
iq>plied  to  this  system.  The  fact  that  this  system  has  two  incident  beams  rather  than  a  single  inci¬ 
dent  beam  (as  is  the  original  Bruning  and  Lo  formulation)  causes  no  difficulties.  On  the  con¬ 
trary,  it  actually  simplifies  the  calculation  by  introducing  symmetry  to  the  problem.  An  efficient 
computer  program  has  been  developed  to  carry  out  these  calculations.  The  complete  details  of 
this  methc^  and  the  computer  program  are  given  in  ref.  [5]. 

23.  NONSPHERICAL  PARTICLES 

The  illustrations  in  Figs,  la  and  lb  show  spherical  particles.  However,  it  is  obvious  that  the 
method  of  images  is  valid  for  nonspherical  particles  as  well.  A  numerical  procedure,  based  on  the 
invariant  imbedding  T-matrlx  method^  has  been  developed  to  calculate  the  scattering  from  the 
resulting  two-particle  system  consisting  of  the  nonspherical  particle  and  its  mirror  image.  At  pre¬ 
sent,  the  computer  program  that  implements  ttese  calculations  is  limited  to  particles  with  dimen¬ 
sions  that  do  not  exceed  about  three  times  the  wavelength  of  the  radiation.  Because  of  this  limi¬ 
tation,  we  did  not  use  this  method  herein. 

2.4  HEMISPHERICAL  PARTICLES 

The  case  of  a  hemlspherically  shaped  particle  is  special.  The  particle  and  the  image  particle 
combine  to  form  a  single  complete  sphere.  The  resulting  scattering  problem  can  be  solved  easily 
and  exactly  by  the  ^>pUcation  of  Mie  theory.  The  only  difference  between  this  problem  and 
problems  ordinarily  encountered  in  Mie  theory  is  the  presence  of  two  incident  waves  rather  than  a 
single  incident  wave.  A  Mie  scattering  program  has  been  modified  to  handle  this  case. 
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3.  FORWARD-SCATTER  MIE  APPROXIMATION 


Young  obtained  good  agreement  between  experiment  and  a  simple  theoretical  model  based  on 
Mie  theory.  Ihe  geometry  of  this  model  Is  illustrated  in  Fig.  2.  An  incoming  beam  of  light  is 
scattaed  by  the  particle.  The  scattering  process  is  assumed  not  to  be  affected  by  the  presence  of 
the  mirror  and  is  calculated  using  ordinary  Mie  theory.  The  light  scattered  by  the  particle  can 
readi  the  detector  by  the  two  paths  shown  in  the  figure.  The  backscatter  ray.  which  proceeds 
directly  from  the  particle  to  the  detector,  is  scattered  through  a  large  angle.  Therefore,  it  is  usu¬ 
ally  much  weaker  than  the  (reflected)  forward-scatter  ray  and  and  can  be  neglected.  Thus,  in  this 
^)|xoximation,  only  the  reflected  forward-scatter  ray  contributes  to  the  signal  received  by  the 
detector.  It  is  flirther  assumed  that  the  presence  of  the  particle  does  not  interfere  with  the  propa¬ 
gation  of  die  reflected  forward-scatter  ray  from  the  mirror  to  the  detector.  Hence,  the  forward 
scatter  Mie  ^iproxlmation  is  simply  a  Mie  theory  calculation  with  the  scattering  angle  measured 
from  the  direction  of  the  specularly  reflected  beam. 


ngute  2.  Geometry  of  Young's  forward-scatter  Mie  theory  model. 
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4.  RESULTS 


Results  obtained  by  exact  theory  and  Youngs  forward-scatter  Mle  theory  are  presented  and  com¬ 
pared.  The  results  are  expressed  as  differential  cross  sections  when  he  scattering  is  from  a  single 
particle  on  a  reflecting  surface  and  as  bidirectional  reflectance  distribution  functions  (BRDF) 
when  the  scattering  is  from  a  distribution  of  particles  on  the  surface.  The  wavelength  in  these  cal¬ 
culations  is  10  pm;  the  incident  light  beam  is  unpolarized  and  directed  normal  to  tlie  reflecting 
surface.  The  differential  aoss  section,  da(6)/dn,  is  a  function  of  the  scattering  angle.  6,  which  is 
measured  with  respect  to  the  normal  to  the  surface.  The  units  of  the  differential  cross  section  are 
micron^/steradian.  Results  are  calculated  for  two  values  of  tlie  particle  index  of  refraction;  n  = 

1.8,  which  is  a  transparent  (nonabsorbing)  material,  and  n  =  2.7  +  O.Slj  (where  j  is  the  imaginary 
unit),  which  is  an  opaque  (absorptive)  material. 

4.1  DIFFERENTIAL  CROSS  SECTIONS 

Figure  3a  shows  the  difTerenlial  cross  secUoii;>  foi  a  small  particle,  diameter  =  1  pm.  The  exact 
cross  sections  are  much  smaller  than  the  Mie  approximation  cross  sections.  The  reason  for  this  is 
easy  to  understand.  The  field  strength  near  the  conducting  surface  is  small  (because  of  tlie  bound¬ 
ary  condition),  so  there  is  no  incident  field  amplitude  available  to  be  scattered  by  the  particle. 

Figure  3b  shows  cross  sections  for  a  large  particle,  diameter  =  50  pm.  The  most  important  point  to 
note  here  is  that  the  angular  structure  of  the  exact  differential  cross  section  is  quite  similar  to  that  of 
the  simple  Mie  theory  differential  cross  section.  The  small-angle  region  is  dominated  by  a  large 
di^action  peak,  which  is  several  orders  of  maguimde  greater  than  the  scattering  values  outside  of 
this  peak.  Trie  angular  width  of  the  diffraction  peak  is  given  aporoximately  by  the  formula  for  the 
first  dark  ring  in  the  diffraction  pattern  of  a  circular  aperture,  6  =  1.22Xyd  (radians),  where  X  is  the 
wavelength  and  d  is  the  particle  diameter.  The  maximum  amplitude  of  the  diffraction  peak  is  about 
a  factor  of  5  greater  than  the  Mie  theory  prediction  for  the  case  shown  in  Fig  3b.  In  the  region 
outside  the  diffraction  peaks,  both  scattering  functions  have  a  qualitatively  similar  behavior 
although  the  exact  solution  has  a  more  pronounced  oscillatory  structure. 

The  most  pronounced  difference  between  the  exact  theory  and  the  Mie  approximation  differential 
cross-section  results  is  evident  in  Figure  4.  Here  the  scattering  angles  are  fixed,  and  the  differential 
cross  sections  are  plotted  vs  the  particle  diameter.  The  exact  results  show  a  large  amplitude  oscilla¬ 
tion  that  does  not  occur  in  the  simple  Mie  theory  results.  The  amplitude  of  the  diffraction  peak 
calculated  by  exact  theory  oscillates  above  and  below  the  difrraction  peak  calculated  by  Mie  the¬ 
ory.  Figure  3b  shows  a  case  where  this  amplitude  is  near  a  local  maximum.  These  oscillations  in 
the  differential  aoss  .sections  could  possibly  have  real  consequences  in  experimental  measurements 
of  BRDF  when  a  lasa  or  other  narrowband  light  source  is  used  to  illuminate  the  sample  contami¬ 
nated  surface.  The  numba  density  of  large  particles  may  not  be  sufficient  to  average  over  the 
effects  of  the  oscillations.  The  amount  of  light  scattered  by  the  large  particles  will  depend  sensi¬ 
tively  on  their  diametas.  This  could  lead  to  large  statistical  fluctuations  in  the  measured  BRDF  for 
different  particle  samples  that  appear  to  have  the  same  size  distribution. 
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Figure  5  shows  similar  plots  lor  the  absorptive  particles.  Here  we  see  that  the  large  oscillations  in 
the  exact  solutions  have  been  damped  out  at  the  smaller  angles  so  there  is  much  better  agreement 
between  exact  theory  and  the  Mle  approximatiou.  However,  as  Figure  5c  shows,  at  larger  scatter¬ 
ing  angles,  serious  differences  still  exist  between  the  exact  and  Mie  results. 

4.2  BRDF 

BRDF  is  the  appropriate  measure  of  scattered  light  from  a  distribution  of  particles  on  a  surface. 
The  distribution  is  characterized  by  the  particle  size  density  distribution  function,  n(x).  This 
function  is  defined  such  that  n(x)dx  is  the  number  of  particles  on  the  surface  per  unit  area  of  sur¬ 
face  that  have  diameters  between  the  values  of  x  and  x  ^  dx.  The  BRDF  of  the  particle-contami¬ 
nated  stuface  is  defined  by  the  relation 


(1) 


This  definition  of  BRDF  does  not  have  the  cos(e)  factor  in  the  denominator.  It  is  essentially  the 
same  definition  used  by  Young*. 

The  calcuIatiOEts  were  carried  out  using  the  particle  size  distribution  that  is  used  to  define  the 
MIL-STD-1246  surface  cleanliness  levels^.  This  distribution  is  defined  by 


LogiN)  =  0.926[LogHx,)-  Log\x)] 


where  N  is  the  number  of  particles  per  ft-  with  diameter  greater  than  x  microns.  The 
quantity  Xi  is  called  the  cleanliness  level.  The  density  function  for  this  distribution  is 


„(;C)  =  1  g52i££W 
X 


This  function  is  negative  for  x  <  1,  which  is  unphysical.  Therefore,  we  arbitrarily  defined  n(x)  = 
0  in  the  range  0  ^  x  <1.  This  causes  no  difficulty  since  these  small  particles  make  negligible 
contribution  to  the  BRDF  for  10  pm  radiation. 
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Before  actually  calculating  the  BRDF  by  evaluating  the  integral  expression  in  Eq.  (1).  it  is  useful 
to  plot  the  Integrand  of  the  BRDF  integral  in  Eq.  (1)  for  selected  values  of  0.  The  integrand  is  a 
density  weighted  cross  section  given  by 


F(x,0)  =  n(x) 


dajx.d) 

dCi 


(4) 


Figure  6  shows  the  BRDF  integrands  for  the  absorptive  particles  for  the  scattering  angles  l^and 
10*.  These  curves  show  the  relative  contributions  of  each  size  particle  to  the  total  BRDF.  At  1° 


Figure  6.  Density-weighted  differential  scattering  cross  section  vs  particle  diam¬ 
eter  at  fixed  scattering  angles;  X  =  10  pm,  n  =  2.7  +  0.5  Ij;  (a)  angle  =  1°,  (b) 
angle  =  10°. 
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this  ftinction  has  a  maximum  at  about  SO  4m,  and  all  the  particles  (but  the  very  smallest )  make 
signiflcant  contributions  to  the  BRDF.  However,  at  10®,  the  integrand  has  a  maximum  at  about  25 
4,m,  and  particle  Sizes  larger  than  60  pm  make  very  little  contribution  to  the  BRDF.  The  relative 
contribution  of  the  various  particle  sizes  to  the  BRDF  is  examined  further  by  calculating  the  upper 
limits  of  the  BRDF  integral  defined  so  that  the  integral  evaluates  to  10,  25,  50,  75,  and  90  percent  of 
the  final  BRDF  value.  These  upper  diameter  limits  are  plotted  as  a  function  of  the  scattering  angle 
in  Figure  7.  The  diagram  shows  that  at  1®,  the  particles  with  diameter  greater  than  100  pm  con¬ 
tribute  half  of  the  BRDF,  but  at  10”  these  same  particles  contribute  almost  nothing. 


Scattering  Angle  (deg.) 


Figure  7.  Upper  particle  diameter  limit  of  BRDF  integral  that  evaluates  to  10, 
25,  50, 75,  and  90  percent  of  the  final  BRDF  value. 
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Figure  8  compares  the  BRDF  values  calculated  by  exact  scattering  theory  with  the  BRDF  values 
calculated  using  the  forward  scatter  ^4ie  theory  and  with  the  exact  BRDF  values  for  hemispherical 
particles.  The  reflecting  surface  h  js  an  assumed  cleanliness  level  of  300  in  these  calculations.  In 
general,  the  exact  BRDF  values  ire  larger  than  the  Mie  approximation  BRDF  values  by  about  a 
factor  of  2.  In  the  angular  range  0®  S  0  <  30°,  Young's  forward-scatter  Mie  approximation  is  seen 
to  actually  be  a  better  approximation  to  the  hemispherical  particles  than  to  the  spherical  particles. 


Figure  8.  BRDF  for  cleanliness  level  300  surface  at  A,  =  10  pm.  Comparison 
of  exact  theory  for  spherical  particles,  exact  theory  for  hemispheres,  and  forward- 
scatter  Mie  theory  for  spherical  particles,  (a)  n  =  1.8,  (b)  n  =  2.7  -t-  0.5  Ij. 
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Figure  9  shows  a  direct  comparison  BRDF  values  for  absorptive  and  nonabsorptive  particles.  As 
one  might  expect  intuitively,  the  nonabsorptive  particles  produce  a  larger  BRDF  than  the  absorp¬ 
tive  particles.  This  seems  reasonable  since  the  absorptive  particles  absorb  part  of  the  energy 
leaving  less  energy  to  be  scattered. 


Hgure  9.  Comparison  of  BRDFs  for  transparent  (n  »  1.8)  and  absorptive 
(n  =  2.7  +  0.51j)  particles. 
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5.  CONCLUDING  REMARKS 


An  exact  theory  of  scattering  from  spherical  particles  on  a  conducting  plane  surface  was  used  to 
calculate  selected  differential  cross  section  and  BRDF  scattering  results.  These  exact  results  were 
compared  witl  results  obtained  by  Young's  forward  scatter  Mie  theory.  The  effect  of  particle 
shape  on  the  BRDF  was  studied  by  comparing  results  for  spherical  particles  with  those  for  hemi¬ 
spherical  particles.  Absorbing  and  nonabsorbing  particles  were  also  compared. 
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